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Abstract 



We calculate the S-wave decay widths and energy shifts for ir + ir~ and tt^K^ atoms 
. in the framework of QCD + QED. The evaluation - valid at next-to-leading order 

in isospin symmetry breaking - is performed within a non-relativistic effective field 
theory. The results are of interest for future hadronic atom experiments. 

o 



Key words: Hadronic atoms, Chiral perturbation theory, Non relativistic effective 
Lagrangians, Isospin symmetry breaking, Electromagnetic corrections 
O ; PACS: 03.65.Ge, 03.65.Nk, ll.10.St, 12.39.Fe, 13.40.Ks 

1 Introduction 

> 

Nearly fifty years ago, Deser et al. [1] derived the formulae for the decay width 
and strong energy shift of pionic hydrogen at leading order in isospin symme- 
try breaking. Similar relations also hold for 7r + 7r~ [2] and tt~K + atoms, which 
decay predominantly into 27r° and n°K°, respectively. These Deser-type rela- 
tions allow to extract the scattering lengths from measurements of the decay 
width and the strong energy shift. The DIRAC collaboration [3] at CERN in- 
tends to measure the lifetime of pionium in its ground state at the 10% level, 
which will allow to extract the scattering length difference [a® — OqI a ^ 5% 
accuracy. The experimental result can then be compared with theoretical pre- 
dictions for the S-wave scattering lengths [4-6] and with the results from other 
experiments [7]. Particularly interesting is the fact that one may determine 
in this manner the nature of the SU(2)xSU(2) spontaneous chiral symmetry 
breaking experimentally [8] . New experiments are proposed for CERN PS and 
J- PARC in Japan [9]. In order to determine the scattering lengths from such 
experiments, the theoretical expressions for the decay width and the strong 
energy shift must be known to an accuracy that matches the experimental 
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precision. For this reason, the ground state decay width of pionium has been 
evaluated at next-to-leading order [10-15] in the isospin symmetry breaking 
parameter 5, where both the fine-structure constant a and (m u — m d ) 2 count 
as 0(6). The aim of the present article is to provide the corresponding for- 
mulae for the S-wave decay widths and strong energy shifts of pionium and 
the ir^K^ atom at next-to-leading order in isospin symmetry breaking. A 
detailed derivation of the results will be provided elsewhere [16]. The strong 
energy shift of the it^K^ atom is proportional to the sum of the isospin even 
and odd S-wave nK scattering lengths + % . This sum [18-22] is sensitive 
to the combination of low-energy constants 2L r 6 + L r 8 [23]. The consequences 
of this observation for the SU(3) xSU(3) quark condensate [24] remain to be 
worked out. 



2 Non-relativistic framework 

The non-relativistic effective Lagrangian framework has proven to be a very 
efficient method to investigate bound state characteristics [12,15,25,26]. The 
non-relativistic Lagrangian is exclusively determined by symmetries, which 
are rotational invariance, parity and time reversal. It provides a systematic 
expansion in powers of the isospin breaking parameter 5. What concerns the 
7T~~K + atom, we count both a and m u — as order 5. The different power 
counting for the 7r + 7r~ and n~K + atoms are due to the fact that in QCD, the 
chiral expansion of the pion mass difference A^ = M 2 + —M% is of second order 
in m u — md, while the kaon mass difference Ak = M 2 + — M 2 starts at first 
order in m u — rrid- In the sector with one or two mesons, the non-relativistic 
ttK Lagrangian is £ NR = C\ + £2- The first term contains the one-pion and 
one-kaon sectors, 

£ 1 = I(E^B 2 ) + fc t^_ A4o + _A_ + _|L + . 

where E = — VA — A, B = V x A and the quantity h = n,K stands for 
the non-relativistic pion and kaon fields. We work in the Coulomb gauge and 
eliminate the A component of the photon field by the use of the equations 
of motion. The covariant derivatives are given by D t h± = d t h± =F ieA h± 
and Dh± = Wh± ± ieAh±, where e denotes the electromagnetic coupling. 
What concerns the one-pion-one-kaon sector, we only list the terms needed to 
evaluate the decay width and the energy shift of the ir~K + atom at order 5 9 ^ 2 
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and <5 4 , respectively, 



C 2 = C[tt[kItt_K + + C 2 (nlKl7r K + h.c) + C 3 ntK$n K + 



(2) 



The ellipsis stands for higher order terms 1 . We work in the center of mass 
system and thus omit terms proportional to the total 3-momentum. The total 
and reduced masses read 



Ei = M„i + M K i, fa 



M w i M K i 
M w i + M K i ' 



+,0. 



(3) 



The coupling constant C[ contains contributions coming from the electromag- 
netic form factors of the pion and kaon, 



(4) 



where (r 2 +) and (r 2 K+ ) denote the charge radii of the charged pion and kaon, 
respectively. The low energy constants C\, . . . , C3 may be determined through 
matching the nK amplitude at threshold for various channels, see section 3. 

To evaluate the energy shift and decay width of the n~K + atom at next-to- 
leading order in isospin symmetry breaking, we make use of resolvents. For a 
detailed discussion of the technique, we refer to Ref. [15]. Here, we simply list 
the results. We use dimensional regularization, to treat both ultraviolet and 
infrared singularities. Up to and including order 5 9 ^ 2 , the decay into tt°K° is 
the only decay channel contributing, and we get for the total S-wave decay 
width 



— — 5- T HokoC 2 
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where k = [2^ (£+ - So - «V+/(2^ 2 ))] 1/2 is of order 5 1 ' 2 . The function 
develops an ultraviolet singularity as d — > 3, 



e„ = A(/i)- 
A(^)=^-3) 



1 + 2 
' 1 



ln^ + ln^± + ^( n ) _^(1) _ I 
n n n 



d-3 



ln47r - T'(l) 



(6) 



1 The basis of operators containing two space derivatives can be chosen such that 
none of them contributes to the energy shift and decay width at next-to-leading 
order in isospin symmetry breaking [16]. 
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with ip{n) = r'(n)/r(n) and the running scale At order <5 4 , the total energy- 
shift may be split into a strong part and an electromagnetic part, according 
to 



AE n = AR h + AE C „ 



(7) 



For the discussion of the electromagnetic energy shift, we refer to section 4. 
The strong S-wave energy shift reads 



a 3 /j, + 
nn 3 



Ci 



2vr 



2, /'0'"0^(2^ 



4tt 2 



+ 0(£ 5 ). 



(8) 



The results for the decay width (5) and energy shift (8) are valid at next-to- 
leading order in isospin symmetry breaking. 



3 Matching the low-energy constants 



The coupling constants Cj can be determined through matching the non- 
relativistic and the relativistic amplitudes at threshold. The coupling C 3 is 
needed at order 5° only. However, we have to determine both C\ and C 2 at 
next-to-leading order in isospin symmetry breaking. The relativistic ampli- 
tudes are related to the non-relativistic ones through 

^ m;lfc (q; p) = 4 h(p)^( P )^(q)^(q)]^ ^r fe (q; p), (9) 

with Ui(p) = (M 2 + p 2 ) 1 / 2 . The 3-momentum p denotes the center of mass 
momentum of the incoming particles, q the one of the outgoing particles. 
The effective Lagrangian in Eqs. (1) and (2), allows us to evaluate the non- 
relativistic n~K + — > tt°K° and n~K + — > n~K + scattering amplitudes at 
threshold at order S. In the isospin symmetry limit, the effective couplings C±, 
C2 and C3 are 

„ 2n / , _\ 2y / 27T 27T , 

61 = — (a^ + a ), C 2 = a , C 3 = — af, (10) 

//+ /i+ //+ 

where the S-wave scattering lengths 2 = 1/3 (a^ 2 + 2a^ 2 ) and Oq = 
1/3(oq 2 — (k) 2 ) are defined in QCD, at m u = m d and M n = M n +, M K = M K +. 
By substituting these relations into the expression for the decay width (5) 
and the strong energy shift (8), one obtains the Deser-type formulae [1,2]. 
We demonstrate the matching at next-to-leading order in S by means of the 
ir~K + — > 7i~K + amplitude. In the presence of virtual photons, we first have 
to subtract the one-photon exchange diagram from the full amplitude, as dis- 
played in Fig. 1. The coupling constant C\ is determined by the truncated 
part T^ 1 , which contains an infrared singular Coulomb phase 6 C as d — > 3, 



dp" and a are normalized as in Ref. [18]. 
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Fig. 1. Non-relativistic it K + — > 7r scattering amplitude. The blob describes 
the vector form factor of the pion and kaon. denotes the truncated amplitude. 

f£ ± (p;p) = e 2i *T* ± ( P ;p), 

At order 5, the remainder T^ 1 is free of infrared singularities at threshold. 
The real part of T^ 1 is given by 

ReT^p; p) = § + 5 2 lnM + 1 ReA^ + 0(p), (12) 
IPI /•*+ 4M 7r +M jr + 

with = Cia7r/i + + o(5), _B 2 = — C\ a/j^/n + o(S) and 



-Re A 



K+ 



±;± 

thr 



= CW1 + 



2tt 



A(/i) - In 



V 2 



2tt 2 



o(5). 



(13) 



Here, the ultraviolet pole term A(/i) is removed by renormalizing the coupling 
C\. The renormalization of C\ eliminates at the same time the ultraviolet 
divergence contained in the expression for the energy shift (8). The calcula- 
tion of the relativistic n~K + — > tt~K + scattering amplitude was performed 
at 0(p 4 ,e 2 p 2 ) in Refs. [20,21]. Both the Coulomb phase and the logarithmic 
singularity in Eq. (12) are absent in the real part of the relativistic ampli- 
tude at this order of accuracy, they first occur at order e 2 p 4 . The quantity 
Re A^fa 1 denotes the constant term occurring in the real part of the truncated 
relativistic threshold amplitude. The coupling constant C2 may be determined 
analogously by matching the non-relativistic ir~K + — > tt°K amplitude to the 
relativistic one at order 8. 



4 Results for the 7r K + atom 



The result for the decay width and strong energy shift are valid at next- 
to-leading order in isospin symmetry breaking, and to all orders in the chiral 
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expansion. We get for the decay width at order <5 9 / 2 , in terms of the relativistic 
it~K + — > 71-°^° threshold amplitude, 



where 



a 3 iJ,lp* n A 2 (l + K n ), A 



TV 



.^JL-Re^F + .(*), (14) 



J + 



K„ 



M„ + A K + M K+ A„ . 
M w + + M K + 

-Aafi + (aQ + a ) 



[a 



o 



1 a 

ip(n) - h In— 

n n. 



+ o(6). 



(15) 



The outgoing relative 3-momentum 



p : = ^x(eim 2 w0 ,m 2 k0 ) 



1/2 



(16) 



with X(x,y,z) = x 2 + y 2 + z 2 — 2xy — 2xz — 2yz, is chosen such that the 
total final state energy corresponds to E n = S + — a 2 fi + /(2n 2 ). The quantity 
Re A^f 1 is calculated as follows. One evaluates the relativistic tt~K + — > n°K° 
amplitude near threshold and removes the divergent Coulomb phase. The real 
part contains singularities ~ l/|p| and ~ ln|p|//i + . The constant term in this 



expansion corresponds to ReA®®^. The normalization is chosen such that 



A = a n + e. 



(17) 



The isospin breaking corrections e have been evaluated at 0(p 4 , e 2 p 2 ) in Refs. 
[21,27]. See also the comments in section 6. 

We now discuss the various energy shift contributions. According to Eq. (7), 
the energy shift at order <5 4 is split into an electromagnetic part AE™ and the 
strong part AE^ in Eq. (8). The electromagnetic energy shift contains both 
pure QED corrections as well as finite size effects due to the charge radii of 
the pion and kaon, contained in A. The pure electromagnetic corrections have 
been evaluated in Ref. [28] for arbitrary angular momentum /. We checked 3 
that the electromagnetic energy shift at order a A indeed amounts to 
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n 4 n 3 (2/ + l) 



+ C(a 5 lna). 



(18) 



3 We thank A. Rusetsky for a very useful communication concerning technical 
aspects of the calculation. 
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Here, the first term is generated by the mass insertions, the second contains 
the finite size effects and the last stems from the one-photon exchange contri- 
bution. The strong S-wave energy shift reads at order <5 4 , 



AE^ 



3, ,2 



±A'(1 + K'), A' 



8tt£. 



-ReA^ + o(5), (19) 



with 



K = -2a/i + (a+ + a 
In the isospin limit, the normalized relativistic amplitude 



1 a 

ipin) - h ln- 

n n 



+ o(6). 



(20) 



A' 



+ a + e', 



(21) 



reduces to the sum of the isospin even and odd scattering lengths. The cor- 
rections e' have been obtained at C(p 4 ,e 2 p 2 ) in Refs. [20,21]. See also the 
comments in section 6. The result for AE± in Eq. (19) agrees with the one 

obtained for the strong energy shift of the ground state in pionic hydrogen 

i±;± 



if we replace fx + with the reduced mass of the n~p atom and Re A thr 
with the constant term in the threshold expansion for the real part of the 
truncated n~p — > n~p amplitude. 

What remains to be added are the vacuum polarization contributions [14,29], 
which are formally of higher order in a, however numerically not negligible. 
The vacuum polarization leads to an energy level shift AE^f c as well as to 
a change in the Coulomb wave function of the tt~K + atom at the origin 
5ip Ktn (0). For the first two energy levels, AE^ C [14,29] is given numerically in 
table 2, section 6. Formally of order a 2l+5 , this contribution is enhanced due 
to its large coefficient containing (fi + /m e ) 2l+2 . The modified Coulomb wave 
function affects both, the decay width and the strong energy shift, see section 
6. 



As discussed in section 6, the electromagnetic contributions (18) are known 
to a high precision. Further, the strong shift in the nP state is very much sup- 
pressed (order a 5 ). A future measurement of the energy splitting between the 
nS and nP states will therefore allow to extract the strong S-wave energy shift 
in Eq. (19), and to determine the combination + <2q of the 7rK scattering 
lengths. The energy splitting between the 2S and 2P states is given by 



AE 2s _ 2p = A£ 2 h + A£ 2 7 - A£ 2 7 + AE™ C - AE™ 

= -1.4±0.1eV. (22) 

The uncertainty displayed is the one in AE\ only. For the numerical values of 
the various energy shift contributions, see table 2 in section 6. 
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5 Results for pionium 



The decay rate and strong energy shift of pionium can be obtained from 
the formulae in Eqs. (5) and (8) through the following substitutions of the 
masses M K + — > M n +, M K o — > M n o and the coupling constants C\ — > c±, 
C 2 -> V2(c 2 - 2c 4 A n ) and C 3 -> 2c 3 [16]. The Cj are the low-energy constants 
defined in Ref. [15]. The S-wave decay width of the tx + ti~ atom reads at order 
5 9 / 2 , in terms of the relativistic 7r + 7r~ — > 7r°7r° threshold amplitude, 



r = — o-V A 2 (I + K "\ 

L 7r,n g^3 "tt^i^tt l x ly w,n) i 
A w = (3q Oq + e,,-, 

= f («o + 2aQ 2 - y (2a° + a 2 ) 

1/2 




1 a 

^(n) - ^(1) h ln- 

n n. 



(23) 



where k = M 2 +/M 2 {) — 1. The quantity .4^ is defined as in Refs. [13,15]. The 
isospin symmetry breaking corrections e n have been evaluated at 0(p A ,p 2 e 2 ) 
in Refs. [13,15,30]. For the decay width of the ground state at order 5 9 ^ 2 , we 
reproduce the result obtained in Refs. [11,13,15]. The electromagnetic energy 
shift AE^ ll is obtained from Eq. (18) through the above mass substitutions 
and A — > 1/3 (r 2 +). Finally, the S-wave energy shift of the 7r + 7r~ atom reads 
at order <5 4 , in terms of the relativistic one-particle irreducible tt + 7t~ — > 7r + 7r~ 
amplitude at threshold, 



n° 



< (i + K tn ) , 



+ 4 



■ip(n) - if>(l) -- + ln- 
n n 



+ o(S), 



(24) 



where A' n is defined analogously to the quantity A' discussed in section 4. The 
isospin symmetry breaking contributions have been calculated at 0(e 2 p 2 ) 
in Refs. [31,32]. For pionium the energy splitting between the 2S and 2P states 
reads 



AK, 2s _ 2p = AE^ 2 + AE™ - + AE™ - AE™ 21 

= -0.59 ± 0.01 eV. 



(25) 



Again the uncertainty displayed is the one in AE\] 2 only. The numerical values 
for the various energy shifts are listed in table 3, section 6. 
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Sh,i 




S h,2 


7r + 7r atom 
ir^K^ atom 


(5.8 ±1.2) • 1(T 2 
(4.0 ±2.2) • 1(T 2 


(6.2 ±1.2) • 1(T 2 
(1.7 ±2.2) • 1(T 2 


(6.1 ± 1.2) • 1(T 2 
(1.5 ±2.2) • 1(T 2 



Table 1 

Next-to-leading order corrections to the Deser-type formulae. 
6 Numerical analysis 

For the S-wave itit scattering lengths, we use the chiral predictions a(] = 0.220± 
0.005 and a 2 = -0.0444 ± 0.0010 [5,6]. The correlation matrix for a° and a 2 , 
is given in Ref. [6]. For the isospin symmetry breaking corrections to the nn 
threshold amplitudes (23) and (24), we use e n = (0.61 ± 0.16) • 10~ 2 and 
e' w = (0.37 ± 0.08) • 10~ 2 as given in Ref. [15] and [32], respectively. For the 
7iK scattering lengths, we use the values from the recent analysis of data 
and Roy-Steiner equations [22], = (0.045 ± 0.012)M~+ and = (0.090 ± 
0.005)M~+. The correlation parameter for a$ and <2q is given in Ref. [22]. 
The isospin breaking corrections to the ttK threshold amplitudes (17) and 
(21) have been worked out in [20,21,27]. Whereas the analytic expressions for 
e and e' obtained in [20,21,27] are not identical, the numerical values agree 
within the uncertainties quoted in [21]. In the following, we use [21] e = (0.1 ± 
0.1) • IQ- 2 M~1 and e' = (0.1 ± 0.3) • lO^AT; 1 . For the charge radii of the pion 
and kaon, we take (r 2 +) = (0.452 ±0.013) fm 2 and (r 2 + ) = (0.363 ±0.072) fm 2 
[33]. 

We obtain for the decay width of the ground state, 

2rv 3 

r 1 = &a?fij,lfa) 2 (l + 6 K , 1 ), r w>1 = —pi t M - a lf (! + > (26) 

where the corrections Sh,i, h = n,K are given in table 1. The strong energy 
shift reads 



A£„ h = -^(4 + ao)(l + 4,n), 
A^;, B = -^(2flg + ag)(l + (27) 

For the first two energy levels, the corrections 5' hn are specified in table 1. 
As mentioned in section 4, these corrections to the Deser-type formulae are 
modified by vacuum polarization, 

5 h , n - 5 Kn + 51% 5' Kn - 5' Kn + 51% (28) 
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ir^K^ 1 atom 




A TP vac r„"\7l 
Ah nl i eS/ \ 


A£ n [eVJ 




n=l, 1=0 


-0.095 


-2.56 


-9.0 ± 1.1 


(3.7 ±0.4) • 10~ 15 


n=2, 1=0 


-0.019 


-0.29 


-1.1 ±0.1 




n=2, 1=1 


-0.006 


-0.02 







Table 2 

Numerical values for the energy shift and the lifetime of the tt^K^ 1 atom. 



tt + tt atom 


A^ m nZ [eV] 




A^, n [eV] 


r w ,„[s] 


n=l, Z=0 


-0.065 


-0.942 


-3.8 ±0.1 


(2.9 ±0.1) • 10~ 15 


n=2, 1=0 


-0.012 


-0.111 


-0.47 ±0.01 




n=2, 1=1 


-0.004 


-0.004 







Table 3 

Numerical values for the energy shift and the lifetime of the ir + Tr~ atom. 



where 

Sh > n ~ "^moT (29) 

Formally, the contribution <5™^ is of order a 2 , but enhanced because of the 
large coefficient containing /j, + /m e . For the ground state, the corrections [14] 
yield 6%% = 0.45 • 10~ 2 and <5™ c = 0.31 • 10~ 2 . The changes in 5^ and S' nA due 
to 5™ are about 5%. For 5' K l however, the correction amounts to 27%. Here, 
we omit the contributions from 5™, because the uncertainties in 5h, n and S' hn 
are much larger than 5™. 



The numerical values for the lifetime T\ = Tf 1 , (r^i = and the energy 
shifts at next-to-leading order in isospin symmetry breaking are given in table 
2 and 3. The energy shifts due to vacuum polarization AE^f c are taken from 
Ref. [14,29]. In the evaluation of the uncertainties, the correlations between 
the S-wave scattering lengths have been taken into account. For the decay 
width and the strong energy shift of the n^K^ atom, the dominant source 
of uncertainty is due to the uncertainties in the scattering lengths and 
do ■ We do not display the error bars for the electromagnetic energy shifts, 
which stem at order a 4 from the uncertainties in (r 2 + ) and (r 2 + ) only. For 
pionium, the uncertainties of AE™ at order a 4 amount to about 0.7%, while 
for the 7r ± ii' =F atom AE™ is known at the 5% level. To estimate the order of 
magnitude of the electromagnetic corrections at higher order, we may compare 
with positronium. Here, the a 5 and a 5 \na corrections [34] amount to about 
2% with respect to the a 4 contributions. 



10 



7 Summary and Conclusions 



We provided the formulae for the energy shifts and decay widths of the 7r + 7r~ 
and tt^K 11 atoms at next-to-leading order in isospin symmetry breaking. To 
confront these predictions with data presents a challenge for future hadronic 
atom experiments. Should it turn out that these predictions are in conflict 
with experiment, one would have to revise our present understanding of the 
low-energy structure of QCD. 
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